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The isospin-dependence of mechanical and chemical instabilities is investi- 



phenomenological equation of state for neutron-rich matter. Respective roles 

of the nuclear mean field and the 2-body stochastic scattering on the evolution 

^ ' of density and isospin fluctuations in either mechanically or chemically unsta- 

C ' 

ble regions of neutron-rich matter are investigated. It is found that the mean 
• <— 1 . 

I field dominates overwhelmingly the fast growth of both fluctuations, while 



the 2-body scattering influences significantly the later growth of the isospin 
fluctuation only. The magnitude of both fluctuations decreases with the in- 
creasing isospin asymmetry because of the larger reduction of the attractive 
isoscalar mean field by the stronger repuslive neutron symmetry potential in 
the more neutron-rich matter. Moreover, it is shown that the isospin frac- 
tionation happens later, but grows faster in the more neutron-rich matter. 
Implications of these results to current experiments exploring properties of 
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neutron-rich matter are discussed. 
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I. INTRODUCTION 



Explosion mechanisms of supernova and related properties of neutron stars are among 
the most interesting topics of modern nuclear astrophysics. Investigations into these ques- 
tions rely critically on the knowledge about the equation of state of isospin asymmetric 
nuclear matter. Thus, one of the most important subjects in nuclear physics is the study of 
the isospin-dependence of the nuclear equation of state using nuclear reactions induced by 
neutron-rich nuclei or radioactive beams . In these reactions transient states of nuclear 
matter with sufficiently high isospin asymmetries as well as large thermal and compressional 
excitations can be created. The rapid progress in recent experiments with rare isotopes has 
therefore made the study of novel properties of extremely isospin-asymmetric nuclear matter 
possible. The planned Rare Isotope Accelerator (RIA) will further enhance the exploration 
of this new frontier dramatically @ . 

Prospects for discovering new physics in neutron-rich matter have generated much in- 
terest in the nuclear science community In particular, the isospin-dependence of the 
nuclear equation of state (EOS), reflected mainly in the density dependence of the symmetry 
energy, has recently received much attention. This is because the isospin-dependence of the 
nuclear EOS is among the most important but very poorly known properties of neutron-rich 
matter [|TU|-|H)[] . It is very important to the mechanisms of Type II supernova explosions and 
neutron-star mergers. It also determines the proton fraction and electron chemical poten- 
tial in neutron stars at (3 equilibrium. These quantities consequently influence the cooling 
rate of protoneutron stars and the possibility of kaon condensation in dense stellar matter 
|16|^0|. Moreover, the isospin-dependence of the nuclear EOS also determines the stability 
boundaries of asymmetric nuclear matter. These boundaries and their isospin-dependence 
subsequently affects the multifragmentation mechanism in nuclear reactions at intermediate 
energies. Multifragmentation is thought to occur in symmetric nuclear matter due to the 
growth of instabilities triggered by density fluctuations, such as the Coulomb, surface and 
volumetric instabilities |H] . While in isospin-asymmetric nuclear matter, it is interesting to 
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study not only the isospin-dependence of the above instabilities, but also new mechanisms 
for nuclear multifragmentation, such as the chemical instability 

In this work we shall explore the isospin-dependence of both mechanical and chemical in- 
stabilities using a Skyrme-type phenomenological EOS for neutron-rich matter within both 
a thermal model and a nuclear transport model. The paper is organized as follows. In 
Section 2, we shall first outline the isospin-dependent nuclear EOS that we use in this work, 
and then investigate relevant thermodynamical properties of neutron-rich matter using a 
thermal model. In Section 3, we shall establish boundaries of the mechanical and chemi- 
cal instabilities and their isospin-dependence. Then, in Section 4 we study the dynamical 
growth of mechanical and chemical instabilities and associated phenomena within a nuclear 
transport model. Finally, a summary and outlook are given in Section 5. 

II. EQUATION OF STATE AND THERMODYNAMICAL PROPERTIES OF 

NEUTRON-RICH MATTER 

To understand the thermodynamic properties of isospin-asymmetric nuclear matter, we 
must have a good knowledge of its EOS. The isospin-dependent part of the nuclear EOS is 
critical for many unique properties of asymmetric matter. At present, nuclear many-body 
theories predict vastly different isospin-dependent nuclear EOS depending on both the calcu- 
lation techniques and the bare two-body and/or three-body interactions employed, see e.g., 
fl~0HT2"|. Therefore, we use here a Skyrme-type phenomenological EOS for asymmetric nu- 
clear matter. With proper parameterizations, it allows us to sample predictions by different 
many-body theories and to calculate many thermodynamical properties analytically. 

Various studies (e.g., [p2] , p3|j ) have shown that the energy per nucleon e(p, 5) in asym- 
metric nuclear matter of density p and isospin asymmetry parameter 

^ = (Pn~ Pp)/{pn + Pp) (1) 

can be approximated very well by a parabolic function in 5. At zero temperature, the e(p, 5) 
can be parameterized as 
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e( P , 5) = ^ + Jt-vT + 3 -eW /3 + S(p) ■ 5\ (2) 

In the above u = p/po is the reduced density, e° F is the Fermi energy and a = —123.6 MeV, 
b = 70.4 MeV and o = 2 corresponding to a stiff nuclear EOS of isospin-symmetric nuclear 
matter. The last term is the symmetry energy whose density dependence is currently rather 



uncertain |10|-4T5|. We adopt here a form of 

S(p) = S (p )-u\ (3) 
which was used by Heiselberg and Hjorth- Jensen in their recent studies of neutron stars 



24| with So(po) = 30 MeV and 7 as a free parameter. We note here that a value of about 



7 = 0.6 was obtained by fitting to the result of variational many-body calculations 24,25 



The 7 parameter is directly related to the K sym parameter of asymmetric nuclear matter 



L sj/m 

via 



K ^ = 9 Po^J Z Ip=po = 9S (po) • 7(7 - 1). (4) 
The values of K sym predicted by many-body theories scatter from about —400 MeV to +466 



MeV(e.g. [22|). Experimental values extracted from studying giant monopole resonances of 
asymmetric nuclei do not constrain the -fT sym parameter either. The reported experimental 
values of K sym are between -566 ± 1350 MeV and 34 ± 159 MeV |B . 



Shown in Fig. 1 are the equations of state with 7 = 0.5 and 1 corresponding to K sym = 
-68 MeV and 0, respectively. The saturation points of asymmetric nuclear matter with 
different 5 are linked with the dashed lines to guide the eye. These two 7 parameters lead 
to rather different saturation points especially for neutron-rich matter. The main features 
of the equation of state obtained with 7 = 0.5 and 7 = 1 are surprisingly similarly to those 
based on the Skyrme Hartree-Fock (SHF) and the Relativistic Mean Field (RMF) models 
respectively. 

It is well known that the symmetry energy has a kinetic and a potential contribution 

S(p) = ^e%ui(2i-l) + V 2 (5) 



where V2 is the potential contribution. The corresponding symmetry potential energy density 



is 



w, 



asy 



V 2P 5 2 



(6) 



and the single-particle potential can be obtained from 



dW, 



asy 



asy 
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±(S u~< - 12.7u 2/3 )S + (51,(7 - 1)« 7 + 4.2m§)(5 2 



(7) 



where "+" and "-" are for q = neutron and q = proton, respectively. Shown in Fig. 2 are 
the symmetry energy and the corresponding symmetry potential with three 7 parameters. 
The magnitude of repulsive (attactive) symmetry potentials for neutrons (protons) increases 
with both the density p and isospin asymmetry 5. The most important effect of v 1 sy is to 
cause the migration of neutrons (protons) from relatively high (low) to low (high) density 
regions leading to the isospin fractionation (distillation) effect, and this effect grows with 
increasing S. 

For asymmetric nuclear matter at a finite temperature T, the nucleon chemical potential 
p q corresponding to the EOS of Eq. |2| is given by p5| 



p q = au + bu a + v q + T 



n 



n=l 



n 



(8) 



where At = 27rh 2 /(m q T) 



1/2 



is the thermal wavelength of a nucleon. The coefficients b n 
are obtained from mathematical inversion of the Fermi distribution function [p8fl . From 
the above chemical potentials for neutrons and protons, the global pressure for asymmetric 
nuclear matter can be obtained from the Gibbs-Duhem relation 



dP 
dp 



dpn 

dp 



(9) 



We separate the result according to 



P Pq ~\~ Pasy ~\~ Pkini 



(10) 



where Pq is the isospin-independent nuclear interaction contribution 



P Q = \a P y + b^-u°+\ (11) 
l a + 1 



Pkm is the kinetic contribution 



Pun = T P {1 + - £ 6 n (^) n [(1 + 5)" +1 + (1 - 5) n+1 ] }, (12) 

/ n=l 4 

and P as3/ is the contribution from the isospin-dependent nuclear interaction Vi 

Pasy = (S olPo u< +1 - 8.5 Po ul)5 2 . (13) 

The total pressure P depends on the isospin asymmetry 6 through both the kinetic pressure 
Pki n and the asymmetry pressure P aS y We compare the relative values of the three partial 
pressures in Fig. 3 for a typical temperature of 5 MeV with 7 = 0.5 and 5 =0.2, 0.4 and 0.6, 
respectively. It is seen that the pressure increases with the increasing isospin asymmetry 
and P aS y has an appreciable contribution to the total pressure. The isothermal spinodal 
(ITS) line moves toward lower densities as the isospin asymmetry increases. Fig. 4 shows in 
more detail the interplay of temperature T, isospin asymmetry 5 and the 7 parameter. In 
each window, the isospin asymmetry 5 goes from to 1 (bottom to top) with an increment 
of 0.2. By increasing the isospin asymmetry the pressure increases as if the temperature 
is increasing. It is also seen that the 7 parameter has a small effect on the total pressure 
especially at high temperature because of the dominating role of the kinetic pressure. 

III. ISOSPIN-DEPENDENCE OF MECHANICAL AND CHEMICAL 

INSTABILITIES 

Based on the EOS obtained from various microscopic theories and phenomenological 



models [p9|-|35|, it has long been predicted that isospin-asymmetric nuclear matter under 



certain conditions can be mechanically or chemically unstable, i.e. 

(dP_ 
[dp 



dP\ 

— — < (mechanical), (14) 



T,6 



or 



85 



) 



P,T 



< (chemical). 



(15) 



In these unstable regions, small fluctuations in p or 5 respectively are expected to grow. 

Mechanical instabilities have been well studied in the past for isospin-symmetric systems, 
here we concentrate on probing its isospin-dependence. Results presented in Fig. 3 and 4 
allow us to explore effects of the temperature, density and 7 parameter on not only the 
pressure itself, but also how these factors influence the mechanical instability. The latter 
happens in the region where the slope of the pressure with respect to density is negative, 
as defined in Eq. [14]. From results shown in Fig. 3 and 4, it is easy to relate the increase 
in temperature with a decreasing probability to observe a mechanical instability. It is seen 
that the mechanical instability region shrinks when either the temperature or the isospin 
asymmetry increases. Above a critical temperature T c determined by the condition 



To be more quantitatively, we show in Fig. 5, the critical temperature as a function of 
the isospin asymmetry 5 with a 7 parameter of 0.5. It is interesting to note that the 
critical temperature decreases with increasing isospin asymmetry. Thus one expects to see 
a shrinking mechanical instability region for increasingly more neutron-rich systems. It is 
important to stress here that the critical temperature T c discussed above is only associated 
with the disappearance of mechanical instabilities. For isospin-asymmetric nuclear matter, 
as we shall discuss in detail in the following, the system can still be chemically unstable 
against density fluctuations and undergo fragmentation above the critical temperature T c . 

To study the chemical instability and its isospin dependence we show in Fig. 6 the 
chemical potential isobars /x(T, P, 5) as a function of 5 for neutrons and protons with a 
parameter 7 = 0.5 at a typical temperature of 5 MeV and pressure of 0.005, 0.02, 0.06, 0.12, 
0.2 and 0.8 MeV/fm 3 , respectively. We note that for corresponding pressures, the isobars 
for neutrons and protons begin from the same chemical potential at 5 = as one expects. 




(16) 



the pressure increases monotonically with density, and the mechanical instability disappears. 



S 



The most interesting feature in this plot is that there exists an envelope of pressures (0.005 
MeV/fm 3 < p < 0.8 MeV/fm 3 ) inside of which chemical instability occurs over a range of 5 
where 

/ dfx n 



\ d5 

or 



< (17) 



>0. (18) 



V 95 , 

Inside this region, the system is unstable against isospin fluctuations. For instance, if sev- 
eral neutrons have migrated into a region of chemical instability due to some statistical 
or dynamical fluctuations in a reaction process, the isospin asymmetry parameter 5 in the 
region will increase and the energy of that region will decrease because of its lower neutron 
chemical potential. To minimize the total energy of the system, it is then favorable for the 
system to have even more neutrons move to the chemically unstable region, thus leading to 
the further growth of the isospin fluctuation. Results at other temperatures show similar 
features. However, the boundary of the instability region shrinks as the temperatures T 
increase, i.e., for increasing T, chemical instability occurs for much less isospin asymmetric 
matter. The temperature dependence and its origin will be more clearly illustrated in the 
following. 

Having discussed separately the mechanical and chemical instabilities, we now compare 
their relative boundaries in the configure space of (T, p, 5). The Gibbs-Duhem relation of Eq. 
H can be used to find boundaries of mechanically unstable regions (ITS: isothermal spinodal) 
in the p — 5 plane for given temperatures. Since the chemical instability condition has to be 
evaluated at constant pressures, the following Maxwellian relation 



8p n \ ( dfi n \ ( dp, n \ ( dP \ 1 / dP " 



(19) 



V 95 ) T,P \ 95 /T,p \ 9 P /T,5 \ 9 P Jt,5 \ 95 Jt, p 

is used to find boundaries of the chemically unstable regions (DS: diffusive spinodal). Shown 
in Fig. 7 are the boundaries of the mechanical (thick lines) and chemical (think lines) in- 
stabilities in the p — 5 plane with 7 = 0.5 at T = 5 (lower window), 10 (middle window) 
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and 15 MeV (upper window), respectively. It is seen that the diffusive spinodal enveloping 
the region of mechanical instability extends further out into the plane; the two regions of 
mechanical and chemical instabilities do no overlap. As the temperature increases, both 
instabilities become less prominent over a more narrow range of densities at smaller isospin 
asymmetries. It is thus possible to observe phenomena due to the chemical instability even 
in reactions using stable nuclei at intermediate energies, which can attain local tempera- 
tures up to 20 MeV and isospin asymmetries up to about 0.4 ||. The above features are 
independent of the parameters of the EOS and are in good agreement with those based on 
more microscopic many-body theories |29|-|3T| . In particular, we found that the variation of 
7 parameter has very little effect on the instability boundaries. Therefore, in the following 
a constant of 7 = 0.5 is used. The results shown in Fig. 7 not only provide the motivations 
but also serve as a guidance for our following transport model simulations. 

IV. EVOLUTION OF DENSITY AND ISOSPIN FLUCTUATIONS IN 
MECHANICALLY OR CHEMICALLY UNSTABLE NEUTRON-RICH MATTER 

Density and isospin fluctuations are expected to grow in the mechanically and/or chem- 
ically unstable regions of asymmetric nuclear matter. Information about the respective 
growth rates of these fluctuations and how they depend on the isospin-asymmtry of nuclear 
matter is currently still rather rare. This information is vitally important for the structure 
and stability of both neutron stars and radioactive nuclei as well as mechanisms of nuclear 
multifragmentation in reactions induced by neutron-rich nuclei. In the following, we study 
the evolution of fluctuations in isospin-asymmetric nuclear matter initialized in the mechan- 
ically or chemically unstable regions within an isospin-dependent transport model IBUU 
]. We note that Baran et al have used a similar approach |33| in studying thermody- 



namical instabilities which can be identified either as a mechanical or chemical instability 



34"] 



The IBUU model uses consistently the isospin-dependent EOS and the corresponding 
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potentials as in the thermal model outlined in Section 2. Moreover, scatterings among 
neutrons and protons are fully isospin-dependent in terms of their total and differential 
cross sections as well as the Pauli blockings. Nucleons are initialized using a Boltzmann 
distribution function in a cubic box of length Lb ox with periodic boundary conditions. The 
box is further divided into cells of 1 fm 3 volume in which the average density p ce u and isospin 
asymmetry 5 ce u are evaluated. We used 10 4 test particles per nucleon in evaluating the p ce u 
and 5 ce u. 

Shown in Fig. 8 is an illustration of the evolution of a system initialized in the mechan- 
ically unstable region with Tj = 5 MeV, 5i = 0.6 and pi = 0.05 fm -3 . For the purpose of 
this illustration a small box of L box =10 fm is used. The most interesting feature shown 
here is the gradually increasing isospin fractionation. This is indicated by the spreading of 
the initial system into regions with p < Pi and 5 > 5i and where p > Pi but 5 < <5j. The 
degree of isospin fractionation can be quantified by using the ratio (N/Z) gas / (N/Z)u qui( i, 
where (N/Z) gas and (N/Z) Hquid is the isospin asymmetry of the low (p/po < 1/8) and high 
(p/po > 1/8) density regions, respectively. Shown in Fig. 9 are the degree of isospin fraction- 
ation as a function of time for a system initialized with T; =5 MeV at a density of pi = 0.05 
fm~ 3 . As seen from Fig. 7, the system is mechanically unstable with an isospin asymmetry 
5 of 0.2 and 0.6, while it is chemically unstable with 5 = 0.9. As one expects from the S 
dependence of the symmetry potential v asy , the isospin fractionation grows faster with the 
increasing S. It is also interesting to note that the isospin fractionation happens later as the 
Si increases. 

The variations of the density p and isospin asymmetry 5 with respect to their initial 
values can be quantified by using, respectively, 



^)-(p 2 -p, 2 ) 1/2 



(20) 



and 



<T S (t) = (5 2 




(21) 
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where the average is over all cells. Shown in Fig. 10 and Fig. 11 are the reduced variation 
in isospin asymmetry ag(t) — 0-5(0) and density <Td(t)/Pi as a function of time for a system 
initialized at Tj = 5 MeV, pi =0.05 fm~ 3 and Si =0.2, 0.6 and 0.9, respectively. As a 
reference, results are also shown for a system initialized with Tj = 15 MeV and Si = 0.9 with 
the dash-dot lines (with this initial conditions the system is both mechanically and chemically 
stable as shown in Fig. 7). As one expects, both the isospin and density fluctuations stay 
almost constant and there is no isospin fractionation at all for this system. While for a 
system initialized in the mechanically (chemically) unstable region with Si =0.2 and 0.6 
(Si =0.9), it is interesting to see that both the isospin and density fluctuations grow faster 
with the decreasing isospin asymmetry Si. This finding is consistent with that found by 
Baran et al in ref. []33|| . Thus, the more neutron- rich system in either the mechanically or 



chemically unstable regions is more stable against the growth of both density and isospin 
fluctuations. This is also responsible for the later start of the isospin fractionation in the 
more neutron-rich systems as shown in Fig. 9. 

Why is the more neutron-rich matter more stable against both the density and isospin 
fluctuations? In the following we try to answer not only this question, but also explore 
seeds of fluctuations and their isospin dependence, and investigate which aspects of nuclear 
dynamics are important in governing the growth of fluctuations. The evolution dynamics 
of asymmetric nuclear matter is governed by the isospin-dependent nuclear mean filed and 
stochastic nuclear scatterings. It is necessary to study their respective roles in generating 
the fluctuations. We start with investigating the role of nucleon-nucleon scatterings. In 
our approach, first-order effects of 2-body stochastic scatterings are included through the 
collision integral of the BUU equation. Second-order effects from the explicitly stochastic 
correction to the collision integral as in the Boltzmann-Langevin model are negelected. 
Thus, t here are two main seeds of fluctuations in our approach, i.e., the initial numerical 
fluctuations from the random sampling of the initial state and the later 2-body stochastic 
nuclear scatterings. Both of them may lead to the growth of fluctuations by propagating 
through the nuclear mean field in the early and later stages of the evolution, respectively. 
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The magnitude of the initial numerical fluctuations in both density and isospin asymmetrty 
is clearly shown in the upper window of Fig. 8. The seeds of fluctuations due to 2-body 
stochastic nuclear scatterings are determined by the isospin-dependent collision dynamics 
in asymmetric nuclear matter. Most importantly, the experimental neutron-neutron cross 
section, as shown in Fig. 12, is only about 1/3 of that for neutron-proton collisions |BB 



and also the final state for neutron-neutron scatterings is more strongly Pauli blocked in the 
more neutron-rich matter. To be more quantitative, we show in Fig. 13 the Pauli blocking 
rate as a function of time. As the Si increases from 0.2 to 0.9 the Pauli blocking rate 
increases from about 35% to 50%. Therefore, there are significantly less nucleon-nucleon 
collisions in more neutron-rich nuclear matter. Shown in Fig. 14 are the average number 
of possible (before using the Pauli blocking) nucleon-nucleon collisions per nucleon N try /A 
and the successful ones (N co u/A) as a function of time. The value of N try /A reflects directly 
the isospin-dependent elementary nucleon-nucleon cross sections shown in Fig. 12. Before 
about 40 fm/c, there is essentially no collision because of the strong Pauli blocking when 
the phase space nonuniformity due to the initial fluctuation is still small. Later, nuclear 
scatterings become important, moreover, they are more frequent with the decreasing isospin 
asymmetry. With Si = 0.2 each nucleon can make about 4 scatterings up to the time of 
100 fm/c, while it can only make about 1 collision with Si = 0.9. Because of the combined 
effects of both the Pauli blocking and the isospin-dependent elementary cross sections, the 
number of nuclear scatterings in nuclear matter with Si = 0.9 is consequently only about 
1/6 of that with Si = 0.2, contributing to the smaller isospin and density fluctuations in 
the more neutron-rich matter. To separate the respective roles of the nuclear mean field 
and the 2-body scatterings, we have performed studies by turnning off the nucleon-nucleon 
collisions in the model. As shown in the upper window of Fig. 15, the collisional seeds of 
fluctuations lead to the significant growth of the isospin fluctuation in the later stage of 
the evolution. However, they have very little effect on the growth of density fluctuations 
as shown in the lower window. A comparison of the results obtained with and without the 
collision integral indicates that the growth of fluctuations are overwhelmingly dominated by 
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the nuclear mean field. This is in agreement with the analysis based on the linear response 
theory by Baran et al |[33| . The observed isospin-dependence of the fluctuations can be 
understood from the interplay between the attractive isoscalar mean field and the repulsive 
symmetry potential for neutrons. The symmetry potential, as shown in Eq. [7], is repulsive 
for neutrons and attractive for protons and their magnitudes increase with the increasing 
isospin asymmetry. Thus, the resultant attractive mean field is weaker in the more neutron- 
rich matter. Because of the small number of scatterings, particularly in the early stage of 
the evolution, the growth of fluctuations is thus mainly determined by the strength and sign 
of the resultant nuclear mean field according to the linear response theory. Based on the 
latter, the magnitude of both fluctuations can grow larger with the increasing strength of 
the attractive resultant mean field, and thus also with the decreasing isospin asymmetry S. 



V. SUMMARY AND OUTLOOK 

In summary, utilizing a phenomenological nuclear equation of state within a thermal 
model, first we explored the isospin dependence of the chemical and mechanical instabilities 
in asymmetric nuclear matter. The diffusive spinodal is found to extend further out into the 
(5 — p) plane than the isothermal spinodal. It is shown that the isospin dependence of the 
nuclear equation of state plays a key role in determining the thermodynamical properties 
of asymmetric nuclear matter. We also investigated the evolution of density and isospin 
fluctuations in mechanically and chemically unstable asymmetric nuclear matter within a 
nuclear transport model. Both the isospin and density fluctuations are found to decrease with 
the increasing isospin asymmetry of the system. Respective roles of the nuclear mean field 
and the 2-body stochastic scattering on the evolution of density and isospin fluctuations in 
either mechanically or chemically unstable regions of neutron-rich matter are investigated. It 
is found that the mean field dominates overwhelmingly the fast growth of both fluctuations, 
while the 2-body scattering influences significantly the later growth of the isospin fluctuation 
only. We have also shown that the isospin fractionation happens later, but grows much faster 
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in the more neutron-rich matter. 

It is well known that more neutron-rich systems are less bound and have smaller satu- 
ration densities as illustrated by the saturation lines in Fig. 1. However, fluctuations and 
their growth are also important for determining the final state of a neutron-rich system, 
such as in the projectile fragmentation in producing exotic beams [37]. Our results above 
indicate that fluctuations actually have a compensating role to the lower binding energy in 
stabilizing neutron-rich system where fluctuations are smaller and do not grow as fast as 
in symmetric ones. Moreover, our results on the isospin fractionation accompanying the 
evolution of fluctuations indicate that the configuration of a more dense, isospin-symmetric 
region surrounded by a more isospin-asymmetric gas as in halo nuclei is a natural result 
of the isospin-dependent nuclear dynamics. Furthermore, we expect that the multifrag- 
mentation and isospin fractionation in nuclear reactions induced by neutron-rich nuclei to 
happen on longer time scales compared to symmetric reactions of the same masses. These 
expectations can be tested by measuring products of multifragmentation. In particular, the 
measurement of neutron-neutron, proton-proton as well as fragment-fragment correlation 
functions, in comparative studies of isospin symmetric and asymmetric nuclear reactions 
will be very useful [BRBII] . 
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FIGURES 




-20 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 
0.5 1 1.5 2 

P/Po 

FIG. 1. The equation of state of isospin-asymmetric nuclear matter with the isospin asymme- 
try 5 spans between and 1 with a 7 parameter of 0.5 (lower window) and 1 (upper window), 
respectively. 
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FIG. 2. Symmetry energy (upper window) and potential (lower window) as a function of re- 
duced density with a 7 parameter of 0.5, 1.0 and 1.5, respectively. The symmetry potential is 
plotted for an isospin asymmetry 5 = 0.2. 
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FIG. 3. Total and partial pressures as a function of density, p, for isospin asymmetry 5 = 0.2, 0.4 
and 0.6 with a temperature T = 5 MeV and 7 = 0.5. The isothermal spinodal line (ITS) is indicated 
with arrows. 
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FIG. 4. Total pressure P as a function of density p for isospin asymmetry 
5 = 0.0, 0.2, 0.4, 0.6, 0.8 and 1.0 (from bottom to top) with a temperature T = 5, 10 and 15 MeV 
and 7 = 0.5, 1.0 and 1.5, respectively. 
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FIG. 5. Critical temperature as a function of isospin asymmetry S. 
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FIG. 6. Chemical potential isobars, for neutrons and protons as a function of isospin asym- 
metry, 5 at temperature, T=5 MeV and density parameter, 7 = 0.5. The isobars are at pressures 
of P = 0.005,0.02,0.06,0.12,0.20 and 0.80 MeV/fm 3 , respectively. 
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FIG. 7. Boundaries of the mechanical (thick lines) and chemical (thin lines) in the den- 



sity-isospin asymmetry plane at a temperature of 5, 10 and 15 MeV, respectively. 
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FIG. 8. An illustration of the evolution of asymmetric nuclear matter in the density-isospin 
asymmetry plane. Each point in the scatter plots represent one cell of 1 fm 3 volume in a cubic box 
of side length L^ ox = 10 fm. 
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FIG. 9. Evolution of the degree of isospin fractionation in a cubic box of side length L^ ox 
fm and density 0.05 fm -3 , temperature T{ = 5 MeV and b% =0.2, 0.6 and 0.9, respectively. 
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FIG. 10. Evolution of the reduced density fluctuation in a cubic box of side length L^ ox = 30 
fm and density 0.05 fm~ 3 . The dash-dot lines are calculated with T{ = 15 MeV and Si = 0.9; while 
the solid lines are calculated with T{ = 5 MeV and Si =0.2, 0.6 and 0.9, respectively. 
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FIG. 12. Experimental neutron- neutron and neutron-proton cross sections as a faction of beam 
energy. 
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FIG. 13. The Pauli blocking rate as a function of time in a cubic box of L^ ox = 30 fm with 
density 0.05 fm -3 , temperature T{ = 5 MeV and b% =0.2, 0.6 and 0.9, respectively. 
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FIG. 14. The average number of possible nucleon-nucleon collisions per nucleon (a) and the 
average number of successful nucleon-nucleon collisions per nucleon (b) as a function of time for 
the same systems as in Fig. 13. 
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FIG. 15. Evolution of isospin (upper window) and density (lower window) fluctuations as a 
function of time. The solid and long dashed lines are results of full calculations, while the dotted 
and dash-dot lines are results obtained by turnning off the 2-body scatterings. 
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